Giant non-adiabatic effects in layer metals: Raman spectra of intercalated graphite explained 
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The occurrence of non-adiabatic effects in the vibrational properties of metals have been predicted since 
the 60's, but hardly confirmed experimentally. We report the first fully ab initio calculations of non-adiabatic 
frequencies of a number of layer and conventional metals. We suggest that non-adiabatic effects can be a feature 
of the vibrational Raman spectra of any bulk metal, and show that they are spectacularly large (up to 30% of 
the phonon frequencies) in the case of layer metals, such as superconducting MgE>2, CaCe and other graphite 
intercalated compounds. We develop a framework capable to estimate the electron momentum-relaxation time 
of a given system, and thus its degree of non-adiabaticity, in terms of the experimentally observed frequencies 
and line widths. 

PACS numbers: 63.20.dk,71.15.-m,63.20.kd,74.70.Ad 
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The adiabatic Born-Oppenheimer approximation (ABOA) 
is the state of the art in first-principles calculations of vibra- 
tional properties in solids. In metals, even if in principle un- 
justified, this approximation generally leads to phonon dis- 
persions in very good agreement with experimental data [ 1 ] . 
Indeed violations of the ABOA are hardly visible in solids. 
Engelsberg and Schrieffer suggested J2] that non-adiabatic 
(NA) effects could lead to a significant renormalization of 
zone-center optical phonon frequencies. An intense effort 
have been devoted to detect measurable NA effects in met- 
als, mainly by Raman spectroscopy ol 0,151], however the de- 
tected NA phonon-frequencies-renormalizations are typically 
smaller than some percent of the adiabatic phonon frequency. 

Recently it has been shown that NA effects are crucial 
to interpret the dependence of Raman spectra on doping 
iS 01 in graphene and on doping (H 0] and on diameter and 
temperature! 10] in nanotubes. However, despite being a cen- 
tral issue for the physics of graphene-based systems, the NA 
phonon frequency shift Ao; (the difference between adiabatic 
and NA phonon frequencies), although measurable, is less 
than 1 % of the adiabatic phonon-frequency. Moreover, these 
systems have lower dimensionality and a peculiar electronic 
structure, it is then unclear whether sizeable NA effects can 
actually be observed in truly three-dimensional bulk metals. 

Layer metallic materials, such as Graphite Intercalated 
Compounds (GICs), are three dimensional (3D) metals pos- 
sessing a considerable anisotropy along one direction. Varia- 
tion of the interlayer distance by intercalation or applied pres- 
sure allows to bridge the gap between two-dimensional mono- 
layer systems (such as graphene) and truly three-dimensional 
systems 11 ill . Thus these metals are ideal to judge the role of 
reduced dimensionality in determining non-adiabatic effects. 
CaC6 is also an 11.5 K superconductor! 12] with an interme 
diate electron-phonon coupling A « 0.83 111 311 . Most interest 



ingly, the recent measurement of CaC6 Raman spectrum [14] 
shows that the phonon frequencies related to in-plane C vi- 
brations are almost 80 cmT 1 larger than those obtained from 
density-functional theory (DFT) adiabatic calculations. This 
result is puzzling since in graphite the adiabatic result for the 
E2 5 phonon frequency is 1577 cm -1 , in excellent agreement 



with the experimental value of 1582 cm -1 . 

In this Letter we develop a first-principles theoretical 
framework to calculate the magnitude of NA effects on zone- 
center optic al-phonons in metals. We identify the general con- 
ditions for having sizeable adiabatic effects in 3D bulk metals 
and analyze the experimental constraints that can hinder the 
observation of NA effects. We demonstrate that the occur- 
rence of NA effects is not limited to reduced dimensionality 
but it is a general properties of metals. We apply our ap- 
proach to GICs, MgB2 and a few bulk metals, finding giant 
NA effects. To the best of our knowledge, our work is the first 
systematic implementation and study of non-adiabatic effects 
in a first principles framework. 

Non adiabatic effects due to the treatment of the electron- 
phonon coupling in the Migdal approximation (neglecting of 
vertex corrections in the adiabatic limit) are usually of the or- 
der sJmjM, where m and M are the electronic and ionic 
mass respectively, and thus generally very small. Engelsberg 
and Schrieffer [2] showed that NA renormalization of adia- 
batic phonon frequencies, a larger effect, unrelated to the ne- 
glecting of vertex corrections, and not of the order of ^Jm/M, 
can be observed if the following two conditions are satisfied: 



|q-vp| <C u> 
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where vp is the Fermi velocity, q is the phonon wavevector, 
u> is the phonon frequency, a = H/t and r is the electron 
momentum-relaxation time (Drude) of the electrons near the 
Fermi surface due to all possible momentum-exchange scat- 
tering mechanisms. 

The first condition is verified for optical phonons of small 
wavevector whose phase velocity is larger than the electronic 
Fermi velocity. Even if in principle Eq. ([1) can be fulfilled in 
any system, provided q is small enough, in practice the pene- 
tration length of the laser light limits the smallest exchanged q 
in Raman experiments. The fact that Eq. (Q]) is experimentally 
difficult to satisfy in 3D metals explains why the observation 
of NA effects in these systems has been somehow disappoint- 
ing. 

The second condition states that the electron momentum- 
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relaxation time must be much larger than the phonon period. 
This implies that the electron band-population dynamic is too 
slow to follow the atomic motion and thus the dynamic is non- 
adiabatic |7[]. Eq. © is usually verified in pure and well crys- 
tallized samples at low temperature. Thus Eq. (Q} is the main 
limitation to the observation of NA effects. 

This limitation can be circumvented considering system 
having small \p along certain directions. For example layer 
metals are usually characterized by a small vp component 
perpendicular to the layers. Since the samples are usually 
cleaved parallel to the layers and Raman experiments are per- 
formed perpendicularly to the freshly cleaved surface, the 
scalar product q- \p is small. It is instead also evident that con- 
dition (HJ is easily verified in 2D or ID systems, since Raman 
experiments are performed with an incident light of wavevec- 
tor q perpendicular to the sample spatial dimension(s) and, 
thus, to \f- In fact, NA effects have been consistently ob- 
served in doped graphene and nanotubes in the last few years. 

NA effects can be taken into account by applying 
time-dependent perturbation theory to DFT Neglecting the 
electron-momentum scattering rate (a = 0), the dynamical 
matrix (£>) for a phonon at q=0 (T) and the first-order varia- 
tion of the electronic charge density (An) are lfl5ll : 
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where n(r) is the charge density, the sum is performed on 
Ak k-points, Z?km,kn = (km|AV sc |kn) is the deformation 
potential proportional to the electron-phonon matrix element, 
AV SC and A 2 V h are, respectively, the first derivative of the 
Kohn-Sham potential and the second derivative of the bare 
(purely ionic) potential with respect to the phonon displace- 
ment. The kernel K(r, r') is the second functional deriva- 
tive of the Hartree exchange and correlation potential respect 
to the densities at r and r'. Finally /k„ is the Fermi func- 
tion for a Bloch state |kn) having band energy e-^ n . The 
NA frequencies has to be computed self-consistenty from 
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) /M. Since in Eqs. © and © only inter- 
band transitions contribute, if hcu N <C |eke — £ko|, '•<?. if the 
phonon energy is much smaller than the direct gap between 
empty (eke) and occupied (eko) states, Vr{w NA ) — T>r(0). 
We verified numerically that in all systems considered in the 
present study this approximation applies and that the results 
obtained with T>r(uj NA ) and T>r(0) are indistinguishable. 

The adiabatic frequency ui A is as usual calculated within 
static perturbation theory, as in Ref. (I]]. Note that uj a ^ 
y/T>r{0)/M since in the adiabatic case the intraband (m=n) 
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FIG. 1: Calculated adiabatic (empty symbols) and non-adiabatic 
(rilled symbols) E2 9 phonon frequencies of several GICs, compared 
to the available experimental data. The full line represents the exper- 
imental values. 



term is present [6] and Aw = 
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uj^ is: 
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■S(e F - e k „) = n(e F )g 2 (e F ), 



where n(e F ) is the density of states at the Fermi 



(5) 
level, 



g 2 (e/) = h\D\ (in ] sin \ 2 /2Muj is the square electron-phonon 
matrix element due to intraband transitions, and g 2 (e F ) is the 
average of g 2 (ef) over the Fermi surface. 

We apply our first principles approach to obtain E2 S non- 
adiabatic phonon frequencies at T for a number of GICs, for 
MgB2 and for bulk hep Ti (see footnote IU6I1 for computa- 
tional details). The results are illustrated in Fig. Q]for GICs 
and are collected in Tab. U for all the systems. In general 
huge Aui values (> 60 cm -1 ) are found in all layer com- 
pounds. The most spectacular non-adiabatic effects are found 
in KC 8 (Auj = 310 cm- 1 , 20% of u A ) and in MgB 2 (Aw = 
230 cm -1 , 30% of uj a ). Even in bulk Ti we find a significant 
shift, Auj = 12 cm -1 , (more than 8% of uj a ) . From Fig. Q] it 
is evident that experimental Raman data 1 11 , 18 , T^, 2CJ Eltl 
in all GICs are in much closer agreement with lu na than with 
ui , whereas in MgB2 Raman data lie in between the two the- 
oretical frequencies. In Ti, contrary to the layer metals, the 
observed Raman frequencies B22I1 are much closer to the adi- 
abatic value. 

To further refine our theoretical model, we now compute the 
phonon frequency in presence of a finite value of the electron- 
momentum relaxation rate a, ui a = uj a + Aw". Auj' 7 can 
be obtained using electron and hole Green functions dressed 
by the interaction (electron-electron, electron-phonon, and 
electron-defect) in the calculation of the phonon self-energy 
and of the dynamical matrix J3I [23[ 0] . Assuming i) a inde- 
pendent of energy and ii) n(e) and <?(e/) constants for bands 
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within ep ±a, Refs. Ji, 23, 24] obtained: 

(hu A ) 2 



(ftur 4 ) 2 +CT 2 



(6) 



Note that when hu A > a, Au 



Aoj and the experimen- 



tally observed phonon frequency is w , i.e. the system is 
in the purely NA regime. On the contrary, for ftw A <C a, 
Alu® = and the measured phonon frequency is ui A , i.e. the 
system is completely adiabatic. 

The linewidth of an optical phonon mode at T is also 
affected by the presence of a finite momentum-relaxation 
rate. The decay of a phonon into non-interacting (undressed) 
electron-hole pairs (a = 0) is forbidden for a zone-center op- 
tical mode if the direct gap is larger than the phonon energy 
{huj <C |eke — £ko|)- This condition being verified in all layer 
systems considered here, thus a zero linewidth should be mea- 
sured [25, 26|. However in all stage- 1 GICs and in MgB<i, 
very large linewidths, hardly explainable in term of anhar- 
monicities, are measured. Analogously to the calculation of 
Lo a , a finite momentum-scattering rate a can be considered in 
the evaluation of the imaginary part of the phonon self-energy 
llH |23, 24, 27]. This leads to the following expression for 
the phonon full linewidth at half maximum due to the phonon 
decay in dressed electron-hole pairs (7^ PC ): 



= fiAw- 
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Eq. © can be used to extract a. Setting Auj a = 
we obtain a from the inversion of Eq. ©. Then a is inserted 
in Eq. (0 to determine 7^ and the results are reported in 
Tab. U The comparison between r )^ PC and the experimental 
Raman linewidths r exp is shown in Fig. [2] The agreement 
is overall very good, except for LiCg and RbCs- Note that 
the experimental linewidth includes all sort of broadening (in- 
cluding inhomogeneous effects), while 7J 5 PC includes only 
the phonon decay into dressed electron-hole pairs. Thus when 
T C xp ~ la FC ^ means that the dominant broadening mech- 



anism is the latter one. This is the case in MgB2 11271 12811 . 
and in most GICs. Moreover, by comparing these results with 
the universal curves of Eqs. © and (0 (Fig. [3}, we deduce 
that our estimate of the momentum-relaxation time is a good 
indicator of the degree of non-adiabaticity. 

As an example, our results quantitatively indicate that 
AIgB2, although being characterized by a huge Aw, has 
a smaller Aw" and a relatively short momentum-relaxation 
time. It is thus a mostly adiabatic system, even though Aw ff is 
still very large. On the contrary, all GICs have relatively small 
ct's, falling in the left part of the universal curves of Fig. [3] and 
are thus mostly non-adiabatic. The result on hep titanium ac- 
tually indicates that all metals might have NA T frequencies 
significantly different from the adiabatic ones, but that they 
cannot be experimentally observed by Raman because the 
condition fl} is not verified in experiments. In other words, 
non-adiabaticity is not a unique property of low-dimensional 
or layer systems, and can occur even in perfectly conventional 







to NA 




Alu 


a 


ry EPC 


r ^ exp 


LiC 6 


1362 


1580 


1595 


218 








70 


UC12 


1492 


1591 


1590 


99 


151 


20 


66 


KC» 


1223 


1534 


1522 


311 


245 


120 


157 


KC24, 


1488 


1571 


1599 


83 








26 


RbC$ (4K) 


1300 


1525 


1490 


225 


558 


163 


100 


RbCs (77K) 


1300 


1525 


1480 


225 


650 


180 


120 




1446 


1529 


1525 


83 


325 


36 


71 


CaCtj (300K) 


1446 


1529 


1511 


83 


761 


68 


111 


SrC 6 


1459 


1530 




71 








BaC 6 


1462 


1521 




59 








MgB 2 (2lK) 


538 


761 


600 


224 


867 


199 


197 


Mg 


122 


123 


122.5 


1 


N/A 


N/A 


N/A 


Ti 


139 


151 


141 


12 


N/A 


N/A 


N/A 



TABLE I: Calculated adiabatic and non-adiabatic E2 9 phonon fre- 
quencies, experimental frequencies, theoretical frequency difference, 
calculated a, calculated EPC linewidth, experimental linewidth of 
several GICs and of MqBi. Experimental data are from Refs. 
EHBEimSlifl] Symbols indicate systems for which exper- 
imental data are available, and are the same as in the figures. For Mg 
and Ti the determination of a from our theory is not applicable. We 
note that since our best NA frequency of graphite is 1578 cm -1 , we 
have upshifted by 4 cm -1 all the adiabatic/NA frequencies of GICs. 
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FIG. 2: Calculated EPC linewidth (Eq. |7]l with respect to the ob- 
served experimental ones. 



metals. Interestingly, despite the structural similarity, in bulk 
Mg no NA effect exists. 

In conclusion, we provide a quantitative first-principles the- 
oretical framework to explain the difference between the re- 
ported experimental E2 9 mode frequencies, and the (ordinary) 
adiabatic calculated ones in several relevant layer metals. We 
have shown that giant non-adiabatic effects occur in these sys- 
tems. Moreover, NA effects are in principle very relevant even 
in bulk metals, although difficult to measure. Finally, we have 
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FIG. 3: Upper panel: HAuia/hAiv (Eq. [6] solid line), compared to 
measured frequencies (symbols). Lower panel: J^ PC /hAuj (Eq.[7] 
solid line), and measured linewidths (FWHM). 



shown that the electron momentum-relaxation time can be ex- 
tracted from Raman peak positions and linewidths and is a 
good indicator of the degree of non-adiabaticity of the system. 

Calculations have been performed at the IDRIS French Na- 
tional Computational Facility under the projet CP9-71387. 
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